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Abstract 

We calculate the orbital magnetization of single and double quantum dots 
coupled both by Coulomb interaction and by electron tunneling. The elec- 
tronic states of the quantum dots are calculated in a tight-binding model and 
the magnetization is discussed in relation to the energy spectrum and to the 
edge and bulk states. We identify effects of chirality of the electronic orbits 
and of the anti-crossing of the energy levels when the magnetic field is varied. 
We also consider the effects of detuning the energy spectra of the quantum 
dots by an external gate potential. We compare our results with the recent 
experiments of Oosterkamp et al. [Phys. Rev. Lett. 80, 4951 (1998)]. 
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I. INTRODUCTION 



The magnetic properties of mesoscopic systems are of recent experimental and theoretical 
interest because they give insight on the electronic structure of the system in a noninvasive 
way, in contrast to the transport studieJil which are strongly dependent on the contacts nec- 
essary to measure the electric current and the voltage. Also, while the transport properties 
of the quantum dots depend on the processes that occur at the Fermi level, the magnetic 
properties are determined by the whole energy spectrum of the mesoscopic system. The 
origin of the orbital magnetism consists in the permanent currents carried by the quantum 
eigenstates in the presence of the magnetic field and of the confinement. The currents are 
related to the energy spectra of the system under consideration. For instance, the chirality of 
the current carried by one electron, which gives the direction of the corresponding magnetic 
moment, can be deduced from the dependence of the corresponding eigenenergy on the mag- 
netic field. In general, the spectral properties of the system determine the sample-dependent 
characteristics of the magnetization. 

The study of orbital magnetism under mesoscopic conditions was initiated 
theoretically,!! especially for predicting the role of the boundaries of the samples. An 
enhancement of the magnetization as compared to the Landau diamagnetism of two- 
dimensional (2D) electrons was demonstrated.!! More recently an exactly solvable model 
of non-interacting electrons has been used for simulating persistent currents and obtaining 
the orbital magnetization. This study has shown that the currents are much more sensi- 
tive to the geometry of the system than the associated magnetization.! Bogacheck et a/.! 
have shown, by analytical calculations, for circular and non-circular quantum dots with 
many (noninteracting) electrons that the oscillations of the magnetization with increasing 
the magnetic field have a hierarchy of three characteristic frequencies, due to the oscillations 
of the energy levels in the vicinity of the Fermi energy, due to the Aharonov-Bohm interfer- 
ence, and due to the de Haas-van Alphen (dHvA) effect, respectively. The importance of the 
electron-electron interaction in quantum dots being generally admitted, the consequences 
for the magnetization of 2D electron systems have been studied for finite and infinite modu- 
lated systems,! and also for noncircular dots.!! In the quantum Hall regime, an enhancement 
due to exchange and correlation was shown both experimentally and theoretically.! Other 
recent measurements of the dHvA oscillations, showing clear sawtooth profiles, have been 
performed by Wiegers et a/.! and by Harris et a/..0 For quantum dots, an ingenuous indirect 
technique was used by Oosterkamp et a/.,0 in order to evaluate the change in the magne- 
tization due to single electron tunneling, from transport measurements. Only very recently 
the magnetization of arrays of quantum dots have been directly measured, but insufficiently 
understood.^ 

Since the magnetic moment of an individual dot is extremely small, the experimental 
endeavor is oriented nowadays towards the study of bigger ensembles in order to measure 
a cumulative effect. In this case, the dots can be coupled to each other either only electro- 
statically, or also exchanging electrons by tunneling, so that, in principle, the magnetization 
of the ensemble of coupled dots may be very different from the scaled magnetization of an 
individual dot. In order to be able to distinguish between these two situations, one has to 
know beforehand the behavior of one single dot and then to identify the coupling effects. 

The aim of this paper is a parallel study of magnetic properties of single and double 
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dots in the presence of the electrostatic couphng, resonant tunnehng, and detuning. The 
double dot is viewed as a coherent quantum-mechanical system separated in two regions 
by a constriction that can be controlled by a parameter. The intra- and inter-dot electron- 
electron interaction is taken into account on the same footing, and the detuning between 
the two regions is realized by applying an external potential (or bias) on only one of them. 
The modifications in the energy spectrum, topology of local currents and the consequences 
for the orbital magnetization will be presented. 



II. THE TIGHT-BINDING MODEL AND ITS RANGE OF VALIDITY 

We describe the 2D electron system in perpendicular magnetic field by a discrete, tight 
binding (TB) Hamiltonian, defined on a rectangular lattice (or plaquette). The lattice 
consists of sites along the x-and M sites along the ?/-direction, separated by an inter-site 
distance a, so that a lattice vector reads r„m = nae^ + maey, with n,m integers. Choosing 
the symmetric gauge for the vector potential, the one-electron spinless Hamiltonian reads: 

H = J2[ enm\n,m){n,m\ + t {e^'^'^'^\n,m){n + l,m| + e"*''"'^|ra, m) (ra, m + 1|) + h.c. . (2.1) 

n,m 

In this expression, |n, m) is a set of ortho normal states localized at the sites r„m and (j) is 
the magnetic flux through the unit cell measured in quantum flux units, (j) — Ba^/(f)Q, (pQ — 
1.43 X 10~^^ Weber. The hopping energy t = h /2'ma? will be considered as the energy unit 
and the lattice constant a will be the length unit. The choice of the boundary conditions 
is essential for the spectral properties of the Hamiltonian ( p.l| ) and the orbital magnetism. 
The cyclic boundaries give rise to the Hofstadter butterfly if the commensurability of the 
geometric and magnetic periods is ensured.lii For the finite system the natural boundary 
conditions are of Dirichlet type, in which case the commensurability condition is not neces- 
sary, and a quasi-Hofstadter spectrum, with edge states filling the gaps, is obtained. Also 
the degeneracies specific to the usual Hofstadter spectrum (at _B 7^ 0) are lifted by the 
presence of the infinite walls.01ii 

Our discrete system can be easily tailored into various shapes, or into several subsystems, 
by removing some inter-site hopping terms, or by imposing infinite barriers at some sites. 
For instance, we can model one single quantum dot, and two or more coupled dots. One 
question is whether a real system can be reasonably described by such a discrete model. 
The interatomic distance for GaAs is a ~ 0.5 nm which means that if we wish to reach the 
atomic resolution for a square dot of linear dimensions L = 100 nm, we need a lattice of 
1000 X 1000 sites. Obviously this requires too much memory and computing power, and 
we must restrict to a smaller number of sites. In reality we need much less sites. For 
instance, a grid containing 20 x 20 sites corresponds to an inter-site distance a ~ 5 nm. In 
addition we have to specify the strength of the magnetic field and the number of electrons. 
Obviously, the magnetic length Ib = {fi/eBy/'^ must be larger than the inter-site distance a, 
or equivalently 

0<l/27r. (2.2) 

Therefore, to describe the square dot of 100 nm x 100 nm by a lattice of 20 x 20 sites, we 
have to consider i? < 20T, which actually covers a wide range of experimental interest. As 
another example, a grid 10 x 10 satisfies the same condition only for smaller fields 5 < 5 T. 
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Since the energy spectrum of the one-band 2D tight-binding model contains N x M 
(nondegenerate) eigenvalues, and is bounded in the interval [— 4t,4t], it can accommodate 
at most N X M spinless electrons. Thus, the condition ( |2.2| ) must be supplemented by 
the requirement kpermi < ^/^ which ensures that the cosine-type tight-binding spectrum 
approximates well the parabola of the quasi-free electrons at low magnetic fields. In terms 
of energies the condition can be written as 

Epermi < 2t , (2.3) 

the Fermi energy being measured from the bottom of the spectrum. Evidently, a finer grid 
means a denser spectrum and a larger number of electrons. The second condition indicates 
that a 20 X 20 plaquette provides a reasonable approximation for a system of about 100 
electrons. In general we shall keep the number of electrons N^, below x M/3. Altogether, 
the above conditions show that the lattice model describes correctly a physical quantum 
dot if the electrons occupy only the states corresponding to the bottom-left corner of the 
quasi-Hofstadter spectrum. 

Several papers used the tight-binding model to describe single quantum dots as we do 
here,00 or groups of coupled quantum dots.llilli (In the later case each dot was associated 
to a single site of the lattice and had therefore no internal structure.) The alternative to 
the discrete model is to consider a quantum dot defined by a continuous confining potential, 
to expand the one-electron wave functions in a set of basis functions and to diagonalize the 
corresponding Hamiltonian matrix. But for the numerical calculations the basis has to be 
truncated to a finite set, which is in fact equivalent to choosing a finite number of sites in 
the tight-binding model. The bigger the basis, the bigger the number of electrons afforded 
in the dot. The natural basis functions are the Laguerre polynomials, and for reasonable 
results one needs at least 2-4 times more basis functions than electrons. On the other hand, 
the computational effort increases exponentially with the size of basis set, such that for a 
non-circular and non-parabolic dot one can hardly go beyond 4-6 electrons. § Instead, in the 
tight-binding model we can consider 50-100 electrons or more, for any shape of the dot. 



III. MAGNETIZATION OF A SINGLE DOT 

In this section we discuss the relation of our tight-binding Hamiltonian ( p.l|) and the 
orbital magnetization of the system. The one-body Hamiltonian can be formally written as 

H= Hnm,n'm'\nm){nm'\ , (3.1) 

nm,n'm' 

where Hnm,n'm' = H*,^i nm ^'^^ ^hc matrix elements of H. We denote its eigenvalues by 
and its eigenstates by The position operator for an electron on the lattice can also be 
written as 

r = ^r„m|nm)(nm| , (3.2) 

nm 

and obviously r„m and \nm) are eigenvalues and eigenvectors of the position operator. Then, 
the operator associated with the current carried by one electron of charge e > is 
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J = -er = —[H, r] = — ^ Hnm,n'm.'{rnm " F^'m' ) I "'"^) (^'"^' I ) (3-3) 



nm,?i'm,' 



and obviously any eigenstate of H has a zero average current, (a|J|a) = 0. The terms of 
the sum in Eq. (|3.3| ) can be interpreted as the currents from the site nm to the site n'm', 



Jnm,n'm' T [ Hnm,n'in' i,^nm Tji'm' ) 1^^) I ~^ Hni^i nmi^n'm' I'nm)!'^ ^^TiTTi] 



(3.4) 



With Eq. (|3.4| ) we can calculate the distribution of the current within the system, while 
Eq. ( p.3| ) is sufficient to define the magnetization operator for one electron, as the standard 
orbital magnetic moment. 



M = ^rx J = ^{xJy- yJ^) e, . (3.5) 



For a given eigenstate a of the Hamiltonian the average magnetization is obtained, by 



combining Eqs. ( p.2|) and ( |3.3| ), as 



2e 

{a\M\a) = —— nm' lm[Hnmn'm'{<^\nrn){n'm'\a)] . (3.6) 

nm,n'm' 

This result does not depend on the way the sites are coupled (nearest neighbors, next- 
nearest neighbors, etc.), nor on the presence of the Coulomb interaction. It holds as long as 
we use a one-body Hamiltonian. In particular, in the Hartree approximation, the interaction 
changes the states \a), but not the current operator Eq. (|3.3|) , which is insensitive to the 
diagonal matrix elements of the Hamiltonian. (However, in the Hartree-Fock approximation 
the interaction becomes visible also in the current operator.) 

At fixed magnetic field, the ground state magnetization Mg is calculated by summing up 
the individual contributions M„ of all occupied eigenstates: 

Mg{B)= J2 (3.7) 

Ec,<Ef 

In the localized representation the current density j(r) reads: 

j(r) = — ^5{v - Vnm) Hnm,n'm'ir - Vn'm')\nm){n'm'\ + h.c. (3.8) 

nm n'm' 

The TB model provides a quick proof of the fact that M^^ defined above coincides with 
(a| {—dH/dB) \a), and thus, from the Feynman-Hellman theorem, with dEa/dcj). Indeed, 
from (lOI ): 

^ = = !!^^[^e'--^|nm)(n + l,m| -ne"'™'^|nm)(n,m + l| ] + h.c. (3.9) 

dB 00 dcp 00 nm 

which gives immediately Eq. ( |3.6| ). The same result can be obtained by defining the mag- 
netization as the the magnetic-moment density (l/2meff)r x j(r) integrated over the whole 
area.0 
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The sign of Ma is determined by the chirahty of the corresponding eigenstate, that is by 
the sign of the derivative dEn/d(j), which can be easily observed from the quasi-Hofstadter 
spectrum. The edge and bulk states coexist at strong magnetic fields, have opposite chi- 
ralities, so that their contributions to the total magnetic moment have different signs, and 
eventually may cancel each other. An example is given in Fig. 1 which shows the magneti- 
zation carried by the first 120 eigenstates of a 20 x 20 plaquette at = 0.1, corresponding 
to a square quantum dot of width about 100 nm, at i? ~ 15 T. 

Since the bulk and the edge states generate intercalated bands and quasi-gaps in the 
spectrum, one may guess an oscillatory behavior of the total magnetization Mg of the quan- 
tum dot as a function of the number of electrons accommodated inside. For an infinite 
system these oscillations become the well known sawtooth profile of the total magnetiza- 
tion, the dHvA oscillations in the regime of strong magnetic fields. Recent experiments have 
=hown them tor extended, theoretically inhnite impL|i but, to our knowledge, nothing 
has been reported yet for finite samples, i. e. when the contribution from the boundaries of 
the system is expected to be important, except theoretical results.i'ii3 

In the low field regime the situation is changed: the spectrum of our system is no longer 
organized into bands and gaps, and the consecutive states may have alternating chiralities. 
Therefore the summations of all contributions must give rise to much faster, eventually 
chaotic oscillations of the magnetization as a function of the number of electrons. This 
case was studied by Shapiro, Hajdu and Gurevichi by calculating the magnetic suscepti- 
bility in the grand-canonical ensemble. Their result - for square geometry - is a rapidly 
oscillating comb-type picture which demonstrates an enhancement compared to the Landau 
=u=ceptMityl Here we calculate M, function of at (p = 0.0001 {B ^ 0.015 T) and 
obtain also a comb-like picture which is shown in Fig. 2. 

From an experimental point of view, it is more interesting to analyze the magnetization 
as a function of the magnetic field B. With increasing B, but for a constant number of 
electrons, the magnetization also oscillates changing the sign when the Fermi level crosses 
regions of different chiralities. The internal mechanism can be understood from Figs. 3(a)- 
(b) which show the quasi-Hofstadter spectrum and the magnetization of a 10 x 10 plaquette 
occupied by 10 electrons. One notices that the spectrum contains numerous anti-crossing 
points where the slope of the energy levels change suddenly, together with the chirality of 
the corresponding current. In fact, at these transitions the states change from edge to bulk 
states, or vice versa. However, for the total magnetization only the anti-crossings at the 
Fermi level are important. The change of sign of the magnetic moment of that state changes 
abruptly the total magnetization. Instead, the anti-crossings below Ep do not change the 
total magnetization because there is always a compensation between two adjacent levels. 
Indeed no compensation occurs at the Fermi level. 

For a large number of electrons, the Fermi level increases to a region with a huge number 
of anti-crossings, yielding fine and superfine oscillations of the magnetization obtained by 
Bogacheck et al io. dot. with 2500 electrorr.1 We approach thi. regime in Figs. 4(c)-(d), for 
100 electrons on a 20 x 20 plaquette. The big oscillations of the magnetization are precursors 
of the dHvA effect occurring at strong magnetic fields, while the small oscillations reflect 
the anticrossing points in the energy spectrum at the Fermi energy. 

Some disorder, if present in the system, would lift the degeneracies but would also 
smoothen the anti-crossing regions. Consequently, the oscillations of the magnetization 



6 



would become also smoother than shown in our Fig. 3, and indeed, the same would hap- 
pen at a finite temperature. We shall see in the next section that the electron-electron 
interactions have a similar effect. 

The analysis of the magnetization based on the spectral properties is done here for a 
square geometry and a tight-binding model. The results obtained in Ref.0 prove however 
that the conclusions are model independent. In that paper, a similar behavior showing 
smooth regions and jumps of the magnetization as a function of the field have been found 
for circular and elliptic dots with 2-5 electrons, in the continuum model based on the Darwin- 
Fock Hamiltonian. It was shown that the positions of the jumps in magnetization can be 
identified also from the maxima of the total energy as a function of magnetic field. 

As function of 1/5, the dHvA oscillations are almost equidistant, with amplitude and 
period depending on the number of electrons N^-, as shown in Fig. 4(a). In the scaled vari- 
ables Mg/Nf. and Nf./B we find that all the curves coincide rather well, which is shown in 
Fig. 4(c). From this scaling property one concludes that the period of the dHvA oscilla- 
tions is A{l/B) = const /A^g. This result differs from the usual problem of dHvA effect in 
metals where the period is proportional to l/Ep, while Ep is considered in good approx- 
imation independent of B. It turns out that for small confined systems this is not true, 
the Fermi energy exhibiting large oscillations with maxima coinciding with the peaks of the 
magnetization, as illustrated in Fig. 4(b). 

The scaling of the magnetization curves suggests that the magnetization per particle of 
a quantum dot can be written as 



where f{x) is an oscillating function in the high field regime and it is nearly constant at low 



In principle the spin may play an important role in the magnetization. However, in the 
present paper we neglect the spin.&il The orbital magnetization in GaAs is enhanced due to 
the effective mass by 14.9, compared to the spin contribution. The spin magnetization is in 
general small at moderate or low magnetic fields, when the Zeeman energy is much smaller 
than the cyclotron energy. In a first approximation the spin contribution is independent 
on the current distribution and consists in relatively weak oscillations around zero, see e. g. 
Fig. 6 of Ref.i 

As mentioned in the Introduction we intended to identify possible effects produced by the 
tunneling and electrostatic coupling between the dots. Such effects are important because 
the magnetization is measurable only for ensembles, and not for individual dots.0 Both 
types of coupling have a considerable influence on the orbital part of the wave functions, 
and consequently on the orbital magnetization. This justifies more attention to the orbital 
magnetization for multi-dot systems. 

IV. DOUBLE DOTS: TUNNELING AND ELECTROSTATIC EFFECTS 

We consider now two coupled quantum dots, and label them 1 and 2. The inter-dot 
tunneling, electron-electron interaction, and detuning introduce specific aspects in the dis- 
tribution of charge and persistent currents which affects the orbital magnetization. The 




(3.10) 



fields. 



7 



double dot is sketched in Fig. 5, and it is considered as a unique coherent quantum system 
described by the Hamiltonian 

H = H^{ct>) + iJ2(0, Vg) + riJi2 + H,i_,i , (4.1) 

where Hi and if 2 correspond to the individual dots and are the same as in Eq. ( |2.1| ) . The two 
dots are coupled both by the tunneling term Hu, and by the electron-electron interaction 
Hel-el which here is considered in the Hartree approximation. In addition we also consider 
a gate potential Vg applied on the second dot, i. e. a detuning parameter, which yields an 
energy offset between the two subsystems, and consists of an extra diagonal term in H2. 
The inter-dot resonant-tunneling term is 

Hl2= J2 Cn^mi,n2m2\nimi) {n2m2\ + h.c. , (4.2) 

nimi ,n2m2 

where the sites (nimi) and {n2m2) belong to the first, and to the second dot, respectively. 
The sites (nimi) and (^277^2) are connected or not if c„jmi,n2m2 = 1 or 0. 

For the Coulomb interaction we use the Hartree approximation. The exchange interaction 
might be important at high magnetic fields (even if the spin is ignored), usually well above 
1 Tesla.i'i As we have checked, in the Hartree-Fock approximation the exchange effects are 
not important for the tight-binding model, in the regime studied in this paper. In a recent 
paper Creffield et al. have studied a continuous square quantum dot with two electrons 
in a weak magnetic field by exact diagonalization, and found a good agreement with the 
tight-binding model in the Hartree approximation.E3 

The Coulomb interaction, in the Hartree approximation, reads 

Hei.ei = UeY. E ^=|nm)(nW|, (4.3) 

nm n'm'{^nm) {jl - Tl'Y + {m - TTl')^ 

where Nnm is the mean occupation number of the site r^^ and has to be calculated self- 
consistently with the energy levels. The Coulomb energy, in units of t, becomes Uc = 
I not K, 1 where we have used the dielectric constant k = 12.4 and the effective mass rries = 
0.067 as for GaAs. However, for Uc = 1 it is technically difficult to obtain the convergence 
of our iterative numerical scheme. Therefore, being in fact interested in qualitative results, 
we use in our calculation a lower value, Uc = 0.4, which still produces a strong perturbation 
of the noninteracting states. If we assume the physical dimension of each quantum dot to 
be 50 nm, and we choose a rectangular plaquette of 10 x 20 sites to model the double-dot 
system, meaning a ~ 5 nm, our Uc corresponds to a Coulomb energy of about 8 meV. 

We first consider no exchange of electrons between the dots, i. e. r = 0, and we focus 
our attention on the effects induced by the electron-electron interaction. In principle, this 
interaction may play a role in the magnetic properties of the dots since it produces rear- 
rangements of the electric charge and of persistent currents. One may distinguish between 
intra-dot and the inter-dot interaction, the latter being the electrostatic coupling between 
dots. We would like to find out which one is more important from the point of view of the 
orbital magnetism. 

It has already been observed that in confined svstems the bulk states are more sensitive 
to the Coulomb interaction than the edge states.ilil By calculating the Hartree spectrum, for 
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a fixed number of electrons, we find that tlie energy distance between tlie bulk-type states, 
which almost form energy 'bands', increases, while the energy of the edge states remains 
almost fixed. Therefore the Coulomb interaction mixes the energies of bulk-type and edge- 
type states, and thus states with opposite chirality become intercalated also in the 'bands'. 
This can be seen in Fig. 6, even in the absence of the inter-dot interaction. A supplementary 
repulsion is added when the inter-dot interaction is taken into account, and further changes 
occur in the spectrum. One also notices that the electron interaction contributes to the 
repulsion of the levels at the anti-crossing points. 

The magnetic moments of the Hartree states are sensitive to these effects and incor- 
porate both the broadening of the bands and intercalation of states with opposite chirality. 
We intend to find out the role played by the electron-electron interaction, in general, and 
by the inter-dot coupling, in particular. Comparing the two curves of Mg versus B in Fig. 6, 
we observe that the inter-dot coupling attenuates the jumps and rounds off the peaks. This 
happens mainly at low fields. The comparison of the magnetization shown in Fig. 3, where 
Uc = 0, and those shown in Fig. 6 indicate that the major interaction effect in the magneti- 
zation comes from the inter-dot coupling. Also, for intermediate and high magnetic fields, 
the saw-tooth shifts slightly, as was already noticed in Ref.i However, the magnetic suscep- 
tibility X = dM/dB\^^Q is not affected by the presence of the electrostatic coupling. We 
shall see below that this is not the case in the presence of tunnel coupling. 

We now consider the role of the tunneling, which means r 7^ 0. To be within the 
tunneling regime, we choose r = 0.4t. The tunneling now lifts the two-fold degeneracy 
of the eigenstates of the double dot, and doublets of states appear now in the spectrum. 
On the other hand, the eigenfunctions and the associated persistent currents penetrate the 
constriction and are distributed spatially over the whole area of the double dot. Under these 
circumstances, the distinction between edge- and bulk-states, which is applicable for single 
dots, becomes improper. Mixed states may also occur as shown in Fig. 7 (middle panel). 

This means that the chirality of a state can no longer be guessed from the localization 
of that state in the bulk or at the edge of the sample. In the absence of the interaction the 
current associated with the lowest-energy eigenstate of the double dot is distributed in the 
middle of each dot, as shown in Fig. 7 (top panel). But the Coulomb repulsion pushes the 
current distribution towards the edges, as shown in Fig. 7 (bottom panel). In the mixed state 
^ 46 the current shows, in different regions, either clockwise or anti-clockwise circulation 
so that the total chirality of the state can only be found by the explicit calculation of the 
corresponding magnetization M46. 

The numerical calculation of the total orbital magnetization, shown in Fig. 8, indicates 
that its dependence on the magnetic field is affected by the resonant coupling only in the 
low-field domain where the sharp peaks are again smeared (for 10 electrons per dot this 
occurs below B ~ 2T). Unlike the case of electrostatic coupling, the resonant tunneling 
modifies the magnetic susceptibility. Indeed, for the situation presented in Fig. 8, the slope 
at i? = remains negative but is smaller in magnitude compared to the r = case. 

V. REDISTRIBUTION OF CHARGE IN DETUNED DOUBLE DOTS 

When one of the dots has a different confinement than the other, or is subjected to a 
supplementary gate potential Vg, its individual energy spectrum is detuned with respect to 
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the spectrum of the other dot. For the sake of definiteness we consider a negative gate 
potential, such that the dot 2 in Fig. 5 gains positive electrostatic energy, —eVg > 0, being 
pushed energetically upwards. We assume each dot is of linear dimension L = 100 nm, and 
is described by a 10 x 10 lattice. 

The role of the detuning was emphasized in connection with the transport properties of 
double quantum dots.B The differences in the energy levels of the two dots can block the 
transfer of electrons. The absence of the resonance condition implies the absence of doublets 
in the spectrum (case discussed in the previous section) and also a strong divergence of 
permanent currents at the constriction between dots. 

By varying continuously Vg, one allows the redistribution of the electrons between the 
two dots. Even for a weak inter-dot tunneling (i. e. r 7^ 0, but small), the magnetization 
behaves in an interesting fashion. In order to understand the physical process, let us consider 
the case of nearly isolated dots. The spectra of the two dots are identical, but shifted by eVg. 
Then by changing the gate potential, a series of resonances occur and at each resonance, 
one electron is transferred from the dot 2 to the dot I. The transfer may be accompanied 
by a change of chirality if, for instance, the electron moves from an edge state to a bulk 
state. Then, at the corresponding value of the gate potential, the orbital magnetization of 
the whole system has a jump, together with the occupation numbers of each individual dot. 
The jump of magnetization can be positive or negative, depending on the chiralities of the 
initial and final states. Both situations are visible in the top panel of Fig. 9. 

At a fixed magnetic field, the energy spectrum of the double dot contains a multitude 
of anti-crossing points between the levels of the dot I, the horizontal lines in the bottom 
panel of Fig. 9, and of the detuned dot 2, the lines with slope about 1. Indeed, in Fig. 9 one 
can see that the jumps of the magnetization, and of the numbers of electrons in each dot, 
occur simultaneously with the condition Ej^ = -Eat+i, Ej^ being the highest occupied level 
of the double dot. The redistribution of the electrons means that the system undergoes a 
transition from the configuration {Ni,N2) to {Ni + 1,N2 — 1), where Ni and N2 are the 
numbers of electrons in the first and the second dot respectively, and iVi + A'"2 = is the 
total number of electrons in the system. Then, the jump of the ground state magnetization 
can be expressed in terms of the magnetic moments carried by the eigenstates of the dots: 

AMg = Mg{Ni, N2) - Mgim + 1 , - 1 ) = - ' (5-1) 

where M^\^+i and M^^^ are the magnetic moments carried by the states |A^i + 1) and |A^2)) 
in the dots 1 and 2, respectively. 

This situation actually occurs only for weak inter-dot coupling. Obviously, a stronger 
coupling spoils the quantization of the number of electrons and of the magnetization, that is 
the steps are less sharp and the plateaus less evident, as shown by the dashed lines in Fig. 9. 

As a function of the magnetic field, the spectrum of a double dot develops a dual aspect 
obtained by the mixing of the two Hofstadter spectra of individual dots with a relative shift 
of eVg. For instance, some eigenvalues originating in the first quasi-gap of one dot overlap 
with the second band corresponding to the spectrum of the second dot. A similar situation 
occurs also for higher energies. The effect of the mixing on the total magnetization is shown 
in Fig. 10 (b) where many secondary peaks can be noticed as compared to Fig. 3. The 
differences are thus produced by the numerous anti-crossing points in the spectrum of the 
double dot which appear at the intersection of the two detuned quasi-Hofstadter spectra. 
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In closing this section we mention that by including the Coulomb interaction we obtained 
similar results, without important qualitative differences. 

VI. CONCLUSIONS 

Our general aim was to demonstrate the correlation between the spectral properties and 
the orbital magnetization of a quantum dot, by considering that the orbital magnetization is 
the sum of the individual orbital magnetic moments of all eigenstates. The electrostatic and 
tunneling coupling of two quantum dots also bring new effects. We have noticed that the 
anti-crossing points in the spectrum have special significance because of the sudden change 
in the chirality of the electronic orbits occurring at these points, with corresponding sign 
changes of the magnetic moments. 

We have shown that interesting conclusions can be deduced even only from the field 
dependence of the Fermi energy. This is because the anti-crossing effect at the Fermi level 
cannot be compensated by the opposite contribution of the adjacent level, which is empty. 
So, we have put forward the argument for which the tunneling transport data can be relevant 
for the magnetic properties. Oosterkamp et a/.0 took advantage of this fact and performed 
an indirect measurement of the changes in the magnetization under resonant conditions in 
double dots. In the experiment, by sweeping the gate voltage Vg at different magnetic fields 
B, one identifies all pair values [Vg, B) when a peak in the tunneling current occurs. Then, 
a change in magnetization should occur, as we have discussed in the previous section. In 
Ref.0 the change 5M is calculated as being proportional to SVg/SB. The sign of 6M can 
be positive or negative as we have also found in Fig. 9 (a). We have stated also that the 
sign is determined bvthe relative chiralities of initial and final states. 

The experimentsliiJ were done in the weak coupling limit when the mixing of the two 
spectra is negligible and, intuitively, one may think in terms of alignment of energy levels. 
In our approach, this condition is not compulsory since the double dot is treated as a coherent 
quantum system. In the case of a larger coupling r the anti-crossings, the depletion of the 
detuned dot and the changes in magnetization become more smooth. 

The electron-electron interaction, mainly the inter-dot component, affects the distribu- 
tion of the persistent currents and attenuates the oscillations of the orbital magnetization 
especially at low fields. Nevertheless, according to our results the magnetic susceptibility is 
not influenced by the electrostatic coupling. A scaling behavior of the quantum dot magne- 
tization as a function of magnetic field at different number of electrons has been identified 
in the dHvA regime. 
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FIGURES 




FIG. 1. Magnetizations for each of the first 120 eigenstates, for a lattice model of 20 x 20 
sites, in the high-field regime, = 0.1 (or ~ 15 T, see text). The negative magnetizations 
grouped in bands correspond to bulk states, the positive ones fill the gaps and correspond to edge 
states. 



13 




FIG. 2. Total orbital magnetization as function of the number of non-interacting electrons in 
the low field regime, (f) = 0.0001 {B ?a 0.015 T). The square dot of dimension L = 100 nm is 
represented in the tight-binding model by a plaquette of 20 x 20 sites. 
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FIG. 3. (a) - The dependence of the spectrum and Fermi energy (with crosses) on the magnetic 
field, for a square 10 x 10 plaquette with 10 non-interacting electrons, (b) - The corresponding 
total orbital magnetization. The jumps in magnetization correspond to the change of the sign of 
dEp/dB. (c) - The energy spectrum for a 20 x 20 plaquette, with 100 non-interacting electrons, 
and (d) - the magnetization. 
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FIG. 4. (a) - The magnetization of a square quantum dot L = 100 nm versus the inverse mag- 
netic field for different numbers of electrons: 100 (curve A), 60 (curve B), 40 (curve C). (b) - The 
Fermi energy for the same numbers of electrons, (c) - The scaled representation : magnetization 
per particle Mg/Ne versus Ne/B\ the three curves a,b,c drops in a single one. 
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FIG. 5. Sketch of a double dot: tunneling may occur through the constriction (the central 
channel) and a gate voltage applied on the second dot produces detuning effects. 
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FIG. 6. (a) - Energy spectrum, in the Hartree approximation, in the presence of intra-dot elec- 
tron-electron interaction, but with no inter-dot interaction. The crosses show the Fermi energy, 
(b) - The same, but with both intra- and inter-dot interaction. The dots are coupled only electro- 
statically, and the number of electrons is iVg = 10 in each dot. (c) - The magnetization for the two 
cases: for intra-dot interaction only, with the dashed line and crosses, and for the total interaction, 
with the solid line. The whole system has 20 x 10 sites, and the interaction parameter is Uc = 0.4. 
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FIG. 7. Persistent currents in double dot. Top: The local current corresponding to the eigen- 
state # 1 at J7c = 0. Middle: The same for the eigenstate # 46. Bottom: The local current 
corresponding to the Hartree eigenstate # 1 at ?7c = 0.4 
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FIG. 8. The total magnetization of the double dot in the presence of the electrostatic coupling : 
with tunnel coupling r = OAt (full line) and without tunnel coupling r = (dashed) . The straight 
lines (1) and (2) are the tangents at the magnetization curves in the low field domain for the two 
cases. 
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FIG. 9. The magnetization (a), number of electrons in the nondetuned dot 1 (b) and a piece 
of the energy spectrum (c) versus the gate potential Vg. The following parameters have been used: 
4> = 0.1 (corresponding to « 4T), r = 0.1 (the full line) and r = 0.4 (the dashed line). The 
Fermi level corresponding to A'" = 18 is also shown. 
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FIG. 10. The lower part of the energy spectrum (a) and the ground state magnetization (b) of 
a detuned double dot at = 5 mV. The Fermi level corresponds to 20 electrons accommodated 
in the system. The tunnel coupling is r = 0.4. 
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